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The geometric factors in the field commutators and spring constants of 
the measurement devices in the famous analysis of the measurability of the 
electromagnetic field by Bohr and Rosenfeld are calculated using a Fourier- 
Bessel method for the evaluation of folding integrals, which enables one to 
obtain the general geometric factors as a Fourier-Bessel series. When the 
space regions over which the factors are defined are spherical, the Fourier- 
Bessel series terms are given by elementary functions, and using the standard 
Fourier-integral method of calculating folding integrals, the geometric factors 
can be evaluated in terms of manageable closed-form expressions. 
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1. Introduction 

The fundamental importance of the famous paper of Bohr and Rosenfeld [|IJ on the mea- 
surability of the electromagnetic field is acknowledged by most physicists but, curiously, the 
paper seems to have been read by only a few at the time of its appearance in the early 
193O's,0 and undoubtedly by still less in the more recent times.! The Bohr-Rosenfeld (BR) 
paper arose in a response to Landau and Peierls ||, who argued that, in principle, the 
electromagnetic field is not measurable in a domain where quantum and relativistic effects 
are important. In their paper, BR have refuted this claim by showing that in quantum 
electrodynamics, just as in nonrelativistic quantum mechanics, there is a complete harmony 
between the theoretical formalism and the physical possibilities of measurement. The essen- 
tial aspects of the BR analysis that enabled them to reach this conclusion are the realization 



*Present address: National Institute for Occupational Safety and Health, 1095 Willowdale Road, 
Morgantown, WV 26505, USA 

1 A. Pais says in his book on N. Bohr M: "From decades of involvement with quantum field 
theory I can testify that nevertheless it has been read by very very few of the aficionados. The 
main reason is, I think, that even by Bohr's standards this paper is very difficult to penetrate. It 
takes inordinate care and patience to follow Bohr's often quite complex gyrations with test bodies... 
As a friend of Bohr and mine once said to me: 'It is a very good paper that one does not have to 
read. You just have to know it exists.' Nevertheless men like Pauli and Heitler did read it with 
great care." 

2 A contributing factor to that must have been that the paper had not been available in an English 
translation for a long time. 
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that only field quantities averaged over finite space-time regions are physically meaningful, 
and the employment of essentially classical test bodies of finite size that can have arbitrarily 
large charge and mass, which removed immediately the limit on the field measurability due 
to the radiation-reaction force on the point test charges employed by Landau and Peierls. 
BR have shown that the electromagnetic-field effects of such finite-size test bodies can be 
minimized to exactly the extent demanded by the formalism's commutation relations by 
using classical spring and lever mechanisms connecting the test bodies to the frame of refer- 
ence and to each other, together with the deployment of neutralizing bodies that occupy the 
same space regions as the test bodies but which are charged oppositely and remain attached 
rigidly to the reference frame during the duration of a measurement. The harmony between 
the possibilities of definition afforded by the quantum-electrodynamic formalism and the 
possibilities of measurement could be attained only by the masterful exploitation by BR 
of all the opportunities offered to measurement by classical physics while at the same time 
paying due regard to the limitations imposed on the latter by quantum mechanics.! 

Several illuminating commentaries on the BR analysis have been written 0J^,^,0, some 
of them by authors who were close to the original controversy. While the BR field mea- 
surement philosophy has not been accepted unreservedly by all the writers |J, the technical 
correctness of the BR analysis has not been disputed. Very recently, however, an analysis 
|| of the BR procedure for the measurement of a single space-time-averaged component of 
the electromagnetic field has drawn a conclusion that no compensating spring mechanism 
is needed in order to measure the averaged field component to arbitrary accuracy when no 
neutralizing body is employed. This work is commented on critically elsewhere ||, using 
the calculational methods developed in the present paper. 

The field commutation relations that BR use as the starting point of their analysis, and 
the spring constants of the mechanisms employed in their measurement procedures, are de- 
fined in terms of geometric factors that are double averages over two finite four- dimensional 
space-time regions. As this amounts formally to an eight-dimensional integration, the cal- 
culation of the value of a field commutator for a finite space-time region, or of a BR spring 
constant, is not simple, even though the dimensionality of the integration is reduced, albeit 
not in a straightforward manner, by the presence of a delta function in the integrand. To 
the present author's knowledge, no calculations of these quantities have been reported in the 
literature yet, apart from those for coinciding spherical space-time regions in the comment 
on reference ||. Clearly, a well-controlled algorithm for the evaluation of the field com- 
mutators for finite space-time regions and the BR spring constants is desirable — not least 
because, as we shall see, these quantities amount essentially to the field effects of extended 
charged bodies, a quantitative assessment of which may be needed in special experimental 



3 In the course of their analysis, BR had to examine the problem of measurement of the basic 
mechanical quantities of position, momentum and energy in more detail than it had been done in the 
previous writings of Bohr and Heisenberg, and established the possibility of repeatable momentum 
and energy measurements that may be of arbitrarily short duration, which were rediscovered some 
30 years later by Aharonov and Bohm H. 



2 



situations.Q 

In the present paper, the Fourier-Bessel method developed for an efficient and accurate 
evaluation of multiple folding integrals [[H]] is adapted for the calculation of general BR 
geometric factors. However, the standard Fourier-integral method for calculating folding 
integrals will turn out to be more advantageous in the special case of spherical space regions, 
as it will enable us to evaluate the BR geometric factors with spherical space regions in closed 
form. 

In section 2 of this paper, the BR field commutators are introduced, and the Fourier- 
Bessel method for calculating the general geometric factors in terms of which the field 
commutators and the BR spring constants are defined is developed. In section 3, geometric 
factors with spherical space regions are considered; these are first calculated using Fourier- 
Bessel expansions, and then evaluated in closed form. In the last section, some calculational 
results are presented and discussed, and concluding remarks are made. 

2. Fourier-Bessel expansions of the BR geometric factors 

The starting point in the BR analysis is the set of electromagnetic-field commutation 
relations for the operators of field averages over finite space-time regions instead of for those 
of the field values at space-time points: 

[£®£W] = mm=*(j&-j&), (i) 



\£P,£^] = mKm=^ 



id 



xy i 



(2) 
(3) 

n\u^\ = -m\£^\ = -*(bm - a. (4) 

Here, £®,HjP, etc. are the electric and magnetic field components averaged over a space- 
time region I of volume V\ and duration Tj, as, for example, 
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and the right-hand-side quantities are purely geometric factors defined in terms of double 
averages over space-time regions I and II: 
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4 For example, direct detection of gravitational waves would require repeated measurements of 
very high accuracy on a single object in a regime where quantum effects are important. In these 
so-called "quantum nondemolition" measurements, experimental precision is pushed to the limits 
set by the principles of quantum mechanics and quantum electrodynamics; under such or similar 
circumstances, the measurement procedures and results of the BR analysis may well be of practical 
relevance [11]. 
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where Vi,n = ViVuTiTu and r is the distance between a space point in the region I and a 
space point in the region II. The remaining commutation relations are obtained from (l)-(4) 
and (6)-(8) by cyclic permutations. 

A BR geometric factor, say C^' 11 ^ , can be written as 

4 UI) = I d*! j T " dt 2 J p 1 (r 1 )dr 1 y p 2 (r 2 )dr 2 U(t,r). (9) 

Here, the time intervals associated with the space-time regions I and II are specified, without 
loss of generality, as (0, Ati) and (T, T+ At 2 ) ,! respectively, while the space regions are given 
by the means of uniform density distributions and p 2 {r 2 ) that vanish outside the space 

regions I and II, respectively, and are each normalized to unit volume; the coordinates r\ 
and r 2 now refer to origins located conveniently inside the regions I and II so that the 
displacement r of a space point of region II from a space point of region I is given as 
r = R + r 2 — t*i, with R the displacement of the origin of region II from that of region I. 
The function U(t, r) in equation (9) is the integrand of the multiple integral defining the 
geometric factor; for the geometric factors (6)-(8), it takes the following forms, respectively: 

= 5(t - r) , (ii) 

ox\oy 2 r 

where units such that the speed of light c = 1 are now used, and t = t 2 — t\. The BR 
geometric factor (9) is cast in the form of a time average of a double-folding integral whose 
integrand involves functions that all have finite space extension, as the densities Pk{fk) 
represent finite regions of space, and the radial range of the function U(t, r) is given by 
r = \r\ — t as it contains the delta function 8(t — r) and its derivatives. As such, the BR 
geometric factors are particularly well suited to evaluation by the means of Fourier-Bessel 
HOI. 



expansions 

To this end, the multipoles Ui m (t,r) of the nonspherical function U(t,r), defined by a 
multipole expansion 

U(t,r) = J £U lm (t,r)i l Y lm (f), (13) 

lm 

are expanded as Fourier-Bessel series in spherical Bessel functions ji(q$r) in a range < 
r < r ex : 

oo 

U lm (t,r) = Y,4Z ) m(<l« ) r), (14) 

n=l 



5 BR use the symbols T\ and Tji for these time intervals, and the symbol At for the durations of 
the momentum measurements at the beginning and end of a field-measurement time interval. 
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where q$r ex are the positive roots of ji(x). The coefficients c^(t) in (14) are given in 
terms of the multipoles Ui m (t,r) as 

cfZ\t) = 4y rU lm (t,r) 3l (q^rydr, (15) 
wk > Jo 



Wn 

where 



^ = "fVMr^f. (16) 
The multipole expansion (13) can thus be written for \r\ < r cx as 

1 oo „ 

U(t,r) = -EE^W / eM^-r)Y lm (q^)d^\ (17) 

^ n lm n =l J 

where q$ are vectors with polar angles and discrete moduli \q$\ = q£\ and where the 
identity 

ji(qr)i l Y lm (r) = —f exp(iq-r)Y lm (q) dq (18) 

is employed. 

In the double-folding integral of equation (9), the value of the function U(t, r) is required 
only when 

\r\ = \R + r 2 - n | <r max = R + R 1 + R 2 , (19) 

where R is the separation of the centres of the two densities, and R\ and R 2 are the radii 
beyond which the uniform densities pi(r*i) and p 2 (r 2 ), respectively, vanish — this is simply 
because the product ^1(^1)^2(^2) in the integrand of (9) is bound to be zero when \r\ > r max . 
Thus, the validity of the Fourier-Bessel expansion (17) in the double-folding integral of 
equation (9) will be guaranteed when the expansion radius r cx > r max . Substituting then 
the expansion (17) with r = R + r 2 — t*i and an expansion radius r cx > r max in equation 
(9), one obtains a Fourier-Bessel expansion for the BR geometric factor C^' 1 ^: 

= JZ E E 4H / Pi(-^ ) )P2(^ ) )exp(i^)- J R)^(^) dq«\ (20) 

47r lm n =l J 



Here, 



n \ 1 /" A 'i rT+M 2 ,, , 



are the time averages of the coefficients (15), and p±(— q$) and p 2 (q$) are the Fourier 
transforms of the densities Pk(vk), 

Pk{q) = J Pk{r k ) exp(iq-r fe ) dr k , (22) 

evaluated at the points q = —q$ and q£\ respectively. 
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The Fourier transforms pk(q) can be expanded in multipoles also, 

PM = PtXM^rnMl (23) 

lk m k 

which are given in terms of the similarly defined multipoles PuL k { r k) °f the densities pk{r k ) 
themselves by 

p£l k (q) = 4^ f°° P%l h (r k )ji k (qr k y k dr k . (24) 

J 

When the multipole expansions (23) are substituted in (20), the products of spherical har- 
monics are expanded in terms of single spherical harmonics, and the identity (18) is used 
agjiin, the Fourier-Bessel expansion of the BR geometric factor C^' 11 ^ can be written finally 

1 oo 

cS' u) = ^E E E^C-iHV^-Mix^)^^??) 

Im l\m,\ n=l 




h k \'\ fx i 

i I pi m 




(q^K)i x Y x ,(R), (25) 



where / = (2Z + 1) 1//2 etc., and the large parenthesis denote 3-j coefficients; p! = mi + m 2 and 
p = m + mi + ^2- In this way, when the multipole expansions of the function U(t, r) and 
the densities pkijfk) are given, the evaluation of the eight-dimensional integral (9) is reduced 
to the evaluation of the one- dimensional integrals (15) for cjj^it) and (24) for p\[% (q$)> 

of the two-dimensional integrals (21) for the time averages c[}^ ', and of the Fourier-Bessel 
expansion (25), where, in principle, only the number of terms in the expansion controls the 
degree of approximation to the exact value of the BR geometric factor Cjj' n \ 

We now turn to the calculation of the Fourier-Bessel coefficients (15) and their time 
averages (21). It turns out that with the functional forms (10)-(12) taken by U(t,r), these 
quantities can be evaluated in terms of elementary functions. Let us first determine the 
multipole expansion of the function UA xx (t,r), given by equation (10). First, we regularize 
the space derivative part of UA xx (t,r) by 

d 2 5(t-r) d 2 5(t-r) 

-hm— — K —^, 26 



where the limit e — > is understood to be taken only after a double integration. This yields 



d 2 5(t-r) J 5"{t-r) 3r + e , 3r + e 




+ / , (t - r ) + / , sA * ~ r ) 



dx\d%2 r e-»o [ r + e r(r + e) 2 r(r + e) 

r) <5(i — r) 
r(r + e) r(r + e) 2 



(x 2 - Xi 



i2 



(27) 



6 Cf. equation (21) of reference [10|, which is reconciled with our equation (25) on noting that our 
displacement r = R+r2 — r\ is defined there as r = R+r\—r2-, and that [pi m {q)\* = (— l) m pi- m (q) 
for a real density p(r). 
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while the time derivatives give 



d 2 5(t - r) 5"{t - r) 



(28) 



Now 



-(x 2 - Xl f = fiY^ir) - y u (f )] 2 = | + y/%Y 2 . 2 (r) - \^Y 2Q {r) + Y 22 (r). (29) 

Using this in (27), then substituting (27) and (28) in (10) and taking the limit e — > 
everywhere except in the S(t — r) term of the monopole component, we obtain the function 
UA xx (t,r) as a multipole sum 



5"(t-r) , e (J(i-r) 
— + hm 



iVoo(f) 




n- 2 (f)- VF2o(r)+y 22 (r-) 



(30) 



The limit e — > was performed in the multipole expansion of UA xx (t,r) everywhere except 
in the S(t — r) term of the monopole component because it can be seen easily that it is only 
with this term that the regularization (26) can contribute to the geometric factor i^ 11 ' . No 
regularization is required in the functions (11) and (12), which turn out not to contain any 
monopole components: 




U A (t,r) = -i 



47T 

UB xy (t,r) = u — 



nt-r) | J'jt-r) | J(t-r) 



5"{t - r) 5'{t - r) 



i^io(f). 



i 2 [Y 2 -2(r)-Y 22 (r)], (31) 

(32) 



The multipoles Ui m (r,t) of the various forms of U(t,r) are found easily from equations 
(30)-(32), and the integrals required for the monopole (1 = 0), dipole (I = 1) and quadrupole 
(I = 2) Fourier-Bessel coefficients (15) are evaluated as follows: 



rr 
JO 



rr 
JO 



S"(t - r) | ^ e 5(t - r) 
r e-»o (r + e) 3 r 



jo(qr)r 2 dr 



-q sin(gt)0(t)0(r cx - t) + A^(q, r cx , t) + ± Jim ^9(t)6(r cx - t), 



5"{t - r) 5'(t - r) 



6"(t-r) | z 6'(t-r) | J(t - r 
+A< 2 >(g,r ex ,t). 



ji(gr)r 2 dr = q cos(qt)Q(t)Q(r cx - t) + A (1) (g, r ex , t), 
j 2 (qr)r 2 dr = gsin(gt)0(t)6(r cx - t) 



(33) 
(34) 

(35) 



Here, 
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A (0) (g, r ex , t) = S(t) - cos(qr cx )5(t - r cx ) - r ex j (qr ex )5'(t - r ex ), (36) 

A (1) (g,r cx ,t) = -sm(qr cx )5(t - r ex ) - r ex ji(gr ex )5'(t - r ex ), (37) 

A (2) (g, r ex , £) = [cos(gr ex ) - j (qr cx ) - j 2 {qr ex ))S(t - r cx ) - r cx j 2 (qr cx )5'(t - r ex ). (38) 

These results are obtained following the rules that govern the use of the delta function and 
its derivatives:@ 

/ f{x)5(x - x ) dx = f(x )Q(x - x 1 )Q(x 2 - x ), (39) 

f{x)5'(x - x ) dx = -f'(x )Q(x - x 1 )Q(x 2 - x ) 

+f(x 2 )S(x 2 - x ) - f(x 1 )5(x 1 - x ), (40) 
f(x)8"(x — xo) dx = f"(xo)Q(xo — xi)Q(x 2 — xo) 

f\x 2 )5{x 2 - x ) + f(x 1 )5(x 1 - x ) + f(x 2 )S'(x 2 - x ) - f(x 1 )S'(x 1 - x ), (41) 



XI 

X2 
XI 



where Q(t) is the Heaviside step function: 0(t) — 1 for t > 0, and Q(t) = for t < 0; and 
utilizing the fact that lim^o ji{ x ) — h m x^o d[xji{x)]/dx = 5qi. 

Using equations (30)-(35), the time averaged Fourier-Bessel coefficients (21) needed in 
(25) for the BR geometric factor A^ 1 ' are thus given by 

4°2n = K O) (sm(gi O) t)e(t)0(r cx -t)) - (A(°)(gi°),r ex ,t)) 

- ^lim(e(t + e )- 3 sin(g(°)t)e(t)e(r cx -t))], (42) 

q-n 

ct:i = -^fln = [tfWffiWMr^-t)) + (A^( q ^,r ex ,t))}; (43) 

for the BR geometric factor i^ 11 ' by 

c£n = ^M^IT) [^ 2) (sm(gi 2) t)e(t)0(r cx -t)) + (A^(q^,r ex ,t))}; (44) 



and for the BR geometric factor B^ 1 ^ by 

h7 ' k 1) (cos(g«t)e(t)e(r cx -t)) + (A( 1 )(g«,r cx ,t))l. (45) 



r (10) -i 




3 



Here, the quantities are given by equation (16), and the angular brackets denote the 
time averaging 

I r Ati rT+At 2 

vw> s zszs/. dh l *>/&-«.)■ («) 



7 As a multidimensional integration with finite integration limits is involved here, care has to be 
taken to use the derivatives of the delta function properly. 
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Let us first evaluate the time averages (sin(gt)0(t)0(r cx — t)) and (cos(gt)0(t)0(r cx — t)). 
A straightforward way of doing so is to integrate by parts in both t 2 an d t±, and to utilize 
the fact that dQ(x)/dx = 5(x): 



/•At! rT+Atz 

/ dfi / dt 2 qsm[q(t 2 -t 1 )]G(t 2 -t 1 )e(r m -t 2 + t 1 ) 
Jo Jt 

= { - cos[q(t 2 - h)]e(t 2 - h)e(r ex - t 2 + h) 

JO L t- 2=T 

+ J t cos[g(t 2 - *i)][(J(* 2 - ti)0(r cx - t 2 + h) - Q{t 2 - ^)5(r ex -t 2 + t^dhjdh 



Ati 
ti=0 



= ^ sin[g(T + At 2 - *i)]6(T + At 2 - *i)0(r ex - T - At 2 + h) 

— / sin[g(T + At 2 - h)]Q(T + At 2 - h)5(r ex — T — At 2 + h) dh 
q Jo 

-- sin[g(T - ^)]0(T - tJQ^ -T + h) *\ 
q *i=° 

+ - / sin[g(T - ti)]0(T - h)5(r eK — T + ti) dh 
q Jo 



+ ^{©(ti - T)0(T + At 2 - h) - cos(gr cx )0(r cx - T + tx)0(T + At 2 - r cx - *i)}dti. 

(47) 

Here, terms with smxS(x) were dropped immediately. This gives for the time average 
(sin(gt)0(t)0(r cx -t)>: 

At x At 2 q (sin(gt)0(t)0(r cx - *)) 

= i{ sin[g(T + At 2 - Ati)]6(T + At 2 - At!)0(r cx - T - At 2 + A*i) 

- sin[g(T + Af 2 )]6(T + At 2 )0(r cx — T — Af 2 ) 

- sin[g(T - A^)]0(T - A£i)0(r ex -T + At,) + sin(gT)0(T)0(r cx - T) 

- sin(gr cx )[0(T + At 2 - r cx )0(At! — T — At 2 + r cx ) - 0(T - r cx )0(At! - T + r cx )]} 

- cos(gr cx )0(At 1 - T + r cx )0(T + At 2 - r cx )[min(At 1 , T + At 2 - r cx ) - max(T - r cx , 0)] 
+0(Ati -T)0(T + At 2 )[min(Ati,T + At 2 ) -max(T,0)]. (48) 

The time average (cos(gt)0(t)0(r cx — £)), calculated in a way similar to that for the time 
average (48), yields 

AfiA* 2 q (cos(gt)0(t)0(r cx - *)) 

= i{ cos[g(T + At 2 - Ati)]6(T + At 2 - At!)0(r cx - T - At 2 + A*i) 

- cos[g(T + At 2 )]0(T + At 2 )0(r cx — T — At 2 ) 

- cos[g(T - Aii)]6(T - A^)©^ — T + A*i) + cos(gT)0(T)0(r cx - T) 

- cos(gr cx )[0(T + At 2 - r cx )0(At! — T — At 2 + r cx ) - 0(T - r cx )0(Atx — T + r cx )] 
+0(T + A* 2 )e(A*i - At 2 — T) — 0(T)0(Atx - T)} + sin(gr cx ) 

x0(Atx - T + r cx )0(T + At 2 - r ex )[min(At 1 , T + At 2 - r cx ) - max(T - r cx , 0)]. (49) 
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The ambiguity that may arise in these expressions when the argument of the step function 
vanishes is removed correctly by the definition O(0) = |. The correctness of the somewhat 
lengthy analytical expressions (48) and (49) was checked by performing the two-dimensional 
integration numerically for cases with the time intervals (0, At,) and (T, T + At 2 ) in all the 
possible logical relations of one to another. 

For the time averages {A^(q, r ex , t)), we need the time averages (S(t— r ex )) and 
(S'(t— r cx )), which are evaluated easily to give 

At,At 2 (5(t - rex)) = 9(A*i - T + r cx )0(T + At 2 - r cx ) 
x [min(Ati, T + At 2 - r cx ) - max(T - r cx , 0)], (50) 
At 1 At 2 (8'(t - r cx )) = 0(T + At 2 - r cx )0(At! — T — At 2 + r cx ) 
-0(T-r cx )0(At 1 -T + r cx ). (51) 

These time averages vanish in the limit r cx — > oo. The e — > term in equation (42) can be 
evaluated by taking the limit e — > already after one time integration, and the result is 

- lim(e(t + e)" 3 sin(gt)0(t)0(r cx - t)) 

= °° {At '7JZ +At2) lM - T +Au - ] - max(T - o)l - (52) 

which equals one half of the g-independent term in (48), which, in turn, equals the average 
(S(t — r cx )) with r cx = of equation (50). 

Using equations (48)-(52) and utilizing the fact that the quantities q$r ex are the roots of 
the spherical Bessel functions ji(x) and thus ji(q$r ex ) = 0, the time averaged Fourier-Bessel 
coefficients (42)-(45) are evaluated finally as follows: 

c (00) 



4:71 



3 AhAhw^q^ 

x { sin[gi 0) (T + At 2 - A*i)]9(T + At 2 - Ah)Q(r ex - T - At 2 + A*i) 

- sin[g(°)(T + At 2 )]0(T + At 2 )0(r cx — T — At 2 ) 

- sin[g(°)(T - At^MT - At^Q^ -T + At,) + sin(g(°)T)0(T)0(r cx - T) 
-igWe(A* 1 -T)G(T + A* 2 )[min(Ai 1 ,T + At 2 ) -max(T,0)]}, (53) 

_(2±2) = _ /3-(20) = ±ig (2±2) _ . /27T I 



2 Axxn Axyn ~ VlSA^Arfgi 2 ) 
x { sin[gi 2) (T + At 2 - Ati)]9(T + At 2 - At!)0(r cx — T — At 2 + Ati) 

- sin[gf(T + At 2 )]0(T + At 2 )0(r cx — T — At 2 ) 

- sin[g( 2 )(T - At!)]e(T - A^)©^ — T + At,) + sin(g( 2 )T)e(T)0(r cx - T) 

- sin(gi 2 V cx ){6(T + At 2 - r cx )0(At! - T - At 2 + r cx ) 

-0(T - r cx )0(Atx — T + r cx ) + 0(A*i - T + r cx )0(T + At 2 - r cx ) 
xr cx 1 [min(At 1 , T + At 2 - r cx ) - max(T - r cx , 0)]} 

+gi 2) 0(Atx -T)0(T + At 2 )[min(At 1 , T + At 2 ) -max(T,0)]}, (54) 
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x{ cos[gW(T + At 2 - Ati)]6(T + At 2 - At!)0(r cx -T-At 2 + At ± ) 

- cos[qW(T + At 2 )]0(T + At 2 )0(r cx — T — At 2 ) 

- cos[ ? «(T - Ai!)]e(T - A^)e(r ex - T + A*0 + cos(g«T)0(T)0(r cx - T) 

- cos(gi 1 V cx )[0(T + At 2 - r cx )0(At! — T — At 2 + r cx ) 
-0(T - r cx )0(Atx — T + r cx )] + 0(T + Ai 2 )6(A*i - At 2 - T) 
-0(T)0(A^-T)}. (55) 

Equations (25) and (53)-(55) furnish a general solution to the problem of finding Fourier- 
Bessel expansions of the representative BR geometric factors (6)-(8) with space regions spec- 
ified by the multipoles (24) of their Fourier transforms. Obviously, the formalism developed 
can easily be used to give Fourier-Bessel expansions of all possible BR geometric factors, 
and not only the representatives (6)-(8), as long as their space regions have well-behaved 
multipole expansions. 

3. BR geometric factors with spherical space regions 

3.1 Fourier-Bessel expansions 

Formula (25) gives a Fourier-Bessel expansion of the BR geometric factor C^' IT> for the gen- 
eral case when the uniform densities Pk(r k ) as well as the function U (t, r) are not spherically 
symmetric functions of space coordinates. Let us assume now that the densities are spheri- 
cally symmetric with radii Rk, Pk(fk) — Pk{^k) = (3/47ri?^)0(i?fc — r^). Such an assumption 
should not entail any serious loss in generality as space regions of practical relevance can 
be approximated by regions of spherical shape, but the main reason for considering spher- 
ical space regions is that it simplifies considerably the formulation of the problem and the 
actual calculations. With spherical space regions, only the l\ — Z 2 = terms contribute in 
the general formula (25), Substituting further p^{q) = (47r) 1 / 2 3ji(gi? / fc)/gi?fc for the Fourier 
transforms of the spherical uniform densities Pkij'k), equation (25) simplifies to 

O oo Jim) 

4 I,n) = E E ^Mq^RiM^^Mq^R^UR)- (56) 

But more importantly, this assumption allows an alternative and simpler formulation based 
on a Fourier-Bessel expansion of one of the spherically symmetric densities, say pi(ri), 
instead of the one based on the Fourier-Bessel expansion of the function U(t, r). 

Let us then expand the uniform density pi{r\) as a Fourier-Bessel series in the spherical 
Bessel functions jo{q n ri) in a range < r\ < r ex : 

oo 

Pi( r i) = E c n3o(q n ri). (57) 

n=l 

Here, q n = nn/r cx , and the coefficients c n are given by 
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Cn = — — / Pi ri jo — ri r dr x = — — 7 j 1 — R,' , 58 

r ex Vr cx / Jo Vr cx / 2r£i?i Vr e 



ex 



where R\ = min(i?i, r cx ), with i?i the density's radius. Since the value of the density pi(ri) 
is needed in the multiple integral (9) that defines the BR geometric factor C^' 1 ^ only when 

n = \R + r 2 -r\ < r lmax = R + R 2 + max(T + At 2 ,0), (59) 

where R is the separation of the centres of the two densities, and R 2 and max(T + At 2 , 0) are 
respectively the radii beyond which the density p 2 (r 2 ) and the function U(t,r) vanish, the 
expansion radius r cx should satisfy the relation r cx > ri max . Substituting then the expansion 
(57) with an expansion radius r ex > r lmax in equation (9), and utilizing the identity (18) for 

3o(q n ri) = ^ J exp(iq n -r 1 )dq„ (60) 

with r 1 = R + r 2 — r, one obtains for the BR geometric factor Cjj' 1 ^ a Fourier-Bessel 
expansion: 

1 00 ( ~ 

Cu' U) = t c np2{q n ) / U{-q n ) exp{iq n -R) dq n , (61) 

47r n=l J 



Here, 



P2(q n ) = 4tt r p 2 (r 2 )j (q n r 2 )rldr 2 = 3 ^^1 (62) 
Jo ~ q n R2 



is the Fourier transform of the uniform density p 2 {r 2 ) with a radius R 2 , and U (q) is the time 
averaged Fourier transform of the function U(t, r): 

I /-At! rT+At 2 r 

U(q) = j o dh J t dt 2 J drU(t,r)exp(iq-r). (63) 

With the multipole expansion (13) of U(t, r) and a further use of the identity (18), the time 
averaged Fourier transform (64) can be written also as a multipole sum 

(j{q)=Y.t J i m {q)i l Y lm {q), (64) 

lm 

where 

47Ti' /" A *l fT+At 2 roo 

U im(q) = -r—^T dtl / dt2 / r 2 drU lm (t,r) 3l (qr). (65) 
At 1 Ar 2 Jo Jt Jo 

Using the multipole expansion (64) and again the identity (18), the Fourier-Bessel expansion 
(61) of Cu' 11 ^ now takes the form of a multipole sum: 

oo _ 

4 I,H) = E E CnP2(qn)Ulm(qn)jl(qnR)Y lm (R). (66) 
lm n=l 



12 



Substituting here for c n and p 2 (<? n ) from equations (58) and (62), respectively, and writing 
q n explicitly as nii/r ex , one obtains finally 

i?/ = min( J R 1 ,r ex ), r cx > R + i? 2 + max(T + At 2 , 0). (67b) 

Note that while different spherical Bessel function roots q$r cx and associated weights 
w$ are required with the different multipolarities / of the function U (t, r) in the formula 
(56) based on the Fourier-Bessel expansion of this function, here only the simple quantities 

= nn/r ex and = rg X /2(n7r) 2 are used for all these multipolarities. For this practical 
reason, when one of the densities is a spherically symmetric function, the formulation based 
on the Fourier-Bessel expansion (57) is preferable to an application of the general formula 
(25) to such a case. Moreover, because the infinite-radius Fourier transform of the function 
U (t, r) is now used, rather than the finite-radius transforms (33)-(35), the expressions for the 
time-averaged multipoles (65) of the Fourier transform will be simpler than those required 
in (42)-(45) for the coefficients cf^ needed in the general formulation. It follows from 
equations (15), (21), (42)-(45), (50)-(52) and (65) that the quantities Ui m (mr /r cx ) needed 
in (67) for the BR geometric factor i^ 11 ' are given by 

Uw X \<l) = (4vr) 3/2 jj[g (sin(^)e(O) - §<*(*)>], (68) 

U%Z\q) = -f-U^iq) = 4^ q (sin(gt)6(t)); (69) 
for the BR geometric factor A& 1 ^ by 

U\±2 (?) = Ti4vr^-g (sin(gt)0(t)>; (70) 
and for the BR geometric factor -B^ 11 ^ by 

(j [ * XV \q) = -47r^g(cos(gt)0(t)). (71) 

Here, the time averages (sin(gt)O(t)) and (cos(qt)Q(t)) are the limits r cx — > oo of the finite- 
radius time averages (48) and (49): 

At 1 At 2 q(sm(qt)e(t)) 

= i{ sin[g(T + At 2 - Ati)]6(T + At 2 - At ± ) - sin[g(T + At 2 )]0(T + At 2 ) 

- sin[g(T - Ati)]6(T - A*i) + sin(gT)0(T)} 

+6(A*i -T)G(T + At 2 )[min(At 1 ,T + At 2 ) -max(T,0)], (72) 
AtiAt 2 g(cos(gt)0(t)) 

= -{ cos[g(T + At 2 - Ati)]0(T + At 2 - Ati) - cos[g(T + At 2 )]0(T + At 2 ) 

- cos[g(T - Atx)]0(T - Ati) + cos(gT)0(T) 

+0(T + At 2 )e(Ah - At 2 — T) — 0(T)0(Atx - T)|, (73) 
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and the quantity {S(t)) in (68) is the time average (50) with r ex = 0, which also happens to 
equal the q- independent term in (72) and twice the e — > term of equation (52). 

The greatly simplified (when compared with those of the general formulation) equations 
(67)-(73) give the Fourier-Bessel expansions of the representative BR geometric factors (6)- 
(8) with spherical space regions; it is gratifying to note that it was possible to express the 
terms of these expansions in terms of elementary functions. In fact, it turns out that the BR 
geometric factors with spherical space regions can be evaluated in closed form, which will 
provide a means of testing the accuracy of the Fourier-Bessel results, but this is most easily 
accomplished starting with the standard Fourier-integral treatment of the multiple folding 
integrals involved. 

3.2 Closed-form evaluation 

Using the standard Fourier-integral method of calculating multiple folding integrals, a BR 
geometric factor Gjj' 11 '* with spherical space regions Pfc(r fc ), 

Cg' U) = j o dh J t dt 2 j Pl {r x )dr x j p 2 (r 2 )dr 2 U(t,r), (74) 

can be written as a multipole expansion, the multipoles of which are Fourier integrals: 

4' D5 = 7 r sE / Px{q)h{q)Ui m {q)3l{qR)q 2 dqY lm {R). (75) 

Here, Ui m (q) are the multipoles (65) of the time-averaged Fourier transform of the function 
U(t,r), given by equations (68)-(73), and 

p k (q) = 4vr f pk(r k )jo(qrk)r 2 k dr k = (76) 
Jo qR k 

are the Fourier transforms of the spherically symmetric uniform densities pk{fk)i which have 
radii R k and are normalized to unit volume. This result follows from the convolution theo- 
rem (see, for example, according to which the Fourier transform of a folding integral, 
like the one in equation (74), is the product of the Fourier transforms of the functions 
appearing in the folding integral; equation (75) is simply the time average of the Fourier 
transformation of the folding integral from the momentum space back to the configuration 
space. In general, Fourier integrals of the type of those in equation (75) have to be evaluated 
numerically, and the great advantage of the Fourier-Bessel formulation is that it replaces 
such numerical integration by the analytical method of a series expansion. However, in the 
case of spherical space regions with the Fourier transforms (76), and with the time-averaged 
Fourier transforms (68)-(73) of the functions U(t,r), the Fourier integrals in (75) can in 
fact be evaluated in closed form, and we present such evaluation in this section. 

The quantities Ui m {q) in (75), as given by equations (68)-(73), are linear combinations 
of terms of the form rj (rq) and a g-independent term for the multipolarities I = 0, 2, and 
of terms of the form Tj-\{rq) and a 1/g term for the multipolarity I = 1, while the Fourier 
transforms (76) of the densities are of the form j\{aq)/aq. We define therefore the integrals 

3h(ax)j h (/3x)j h (7x)j k (5x) — , (77) 



ji 4 (n; l 1 ,l 2 ,h,h;a,P,'y,5) 
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the evaluation of which is required in (75) for the parameter values (i) n = 0, l\ = li = 1, 
^3 = 0, U = 0, 2; (ii) n = 0, li = 1% = 1, I3 = — 1, 1 4 = 1; and, on account of the 1/g term in 
(73), (hi) n = 1, h = l 2 = 1, h = (with 7 = 0), Z 4 = 1. 

In principle, it should be possible to evaluate the integral ji 4 (n; Z 1; Z 2 , ^3, ^4; a, /?, 7, 5) in 
closed form for all the integer values n and h,h,h,U for which the integral exists, as the 
integrand can be written as a sum of a finite number of terms, with each term having 
the form asm(hx)/x k or acos(bx)/x h , where k is an integer, and the indefinite integrals 
of such terms can be done in terms of the sine or cosine integrals. In practice, however, 
such a procedure is prohibitively lengthy for all but very small values of the parameters 
n, li, l 2 , 13, 14. for example, the integrand of the integral ji 4 (0; 1, 1, 0, 2; a, (3, 7, 5) has already 
96 terms of the above mentioned form with k ranging from 4 to 8, and the indefinite integral 
of each of these terms generates in turn k new terms, giving in total 572 terms. Furthermore, 
complications arise when the limit x — > (or x — > 00, depending on the definition of the sine 
integral) is taken in the sine integrals, which have arguments of the form ax, as the signs of 
the parameters a, which are linear combinations of the parameters a, (3, 7 and 5, must be 
determined suitably. Remarkably, the computing system Mathematica |13| is able to perform 
the definite integrations in the integrals ji 4 (n; li, 1%, I3, Z4; a, (3, 7, S) required directly. After 
considerable simplifications, and using an economic way of writing linear combinations that 
contain terms with permuting signs by means of the definitions 

a n = a, & =(-!)"& 7n = (-l)K2] 7; ^ = (-1)™% 5' n = (-1)^6, (78) 
where [x] is the integer part of x, the results are as follows: 

u 4 0; 1, 1, 0, 0; a, (3, 7, 5) = ^ 

x [Aa 2 n + {A(3 n - T „ - 5 n ){-3a n + (3 n + ln + <J n )]; (79) 

ji 4 (0; 1, 1, 0, 2; a, (3, 7, 5) = - „ ggn „ x2 2^ 7, j? 

26 880a/3d 2 ^ a n f3 n >y n d n 

x{(a n + (3 n + 7n)(«n + /3„ + 7n - 35„)[6a^ + (6/3 n - 7„)(-5a„ + (3 n + 7n)] 

+ [8(a^ _ I2a n /5 n + (3 2 n ) + 9(a n + (3 n ) ln + >fiK + [24(a n + (3 n ) - ll 7 n]^ - 8#}; (80) 

■•{nil 11 a x\ 71 |«n + /^n + 7n + ^n| 3 

ji 4 (0 1,1,-1,1 a, /3, 7, d) = - jz^— 2_, -r-r 

11 520a/37() ^ a n /* n d„ 

x {5a 3 - 3a^(5/3 n - 37„ + 55 n ) + (-3a n + (3 n + 7„ + <*„) 

x M + (7n - 55 n )(A(3 n - 7 n - *„)]}. (81) 

Here, it is assumed that the parameters a, f3, 7 and 8 are all nonzero; the required integrals 
that have some of these parameters equal to zero were evaluated separately: 

ji 4 (0; 1, 1, 0, 0; a, (3, 0, 0) = -f- £ K+ / ral (a^ - a n /5 n + /£); (82) 

12a/? ^ a n (3 n 
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• /a -i -. a a a n\ 71 \^ ( a n + Pn i In 

ji 4 (0; 1, 1, 0, 0; a, (3, 7, 0) = 2^ 



x[3K-/5 n ) 2 + 2(a n + /5 n )7n-7']; (83) 



ji 4 (0; 1,1,0, 2; a, /3, 0, 5) = - ^ 



x(a n + /3 n - 2 (c4 - 4a„/5„ + /3 2 - S' n 2 ); (84) 

|a n + /3 n + <5'| 3 



ji 4 (l;l,l,0,l;a,A0, ( 5) = - V 

1152a/*) n=0 tt„p n c>; 

x [q£ - 3a„(^ - 4/3^ + 5' n 2 ) + + 5' n )(-3a 2 n + (3 2 n - A(3j n + 5' 2 )]. (85) 

Here, the definition jo(0) = 1 is assumed in the integrals (77); as \im x ^ ji(x) = 
when / > 0, the integrals ji 4 (n; 1, 1, l 3 , U; a, (3, 7, 0) with Z 4 > vanish. The integral 
ji 4 (0; 1, 1, 0, 0; a, (3, 0, 5) that is required also is given already by the integral (83): 

ji 4 (0; 1, 1, 0, 0; a, (3, 0, 5) = ji 4 (0; 1, 1, 0, 0; a, (3, 5, 0). (86) 

Using equations (68)-(73), (75) and (76), the integrals (79)-(86), and the definitions 

a^ ] = W 2 \ 4lT } = 0, a^ = -^-a^=A^ (87) 
4^ = 41^ = 0, = =Fi4^g, b[^ = -A^- (88) 

r = 0, n = T + At 2 - A*!, r 2 = T + At 2 , r 3 = T, r 4 = T-At i; (89) 

^ = "^2A^A^ [min(Atl ' T2) " max ( T3 ' ^ ^ 9 °) 
^ = aT^I ^)©^) " ©(r 3 )e(-r 4 )], (91) 

# = 9( ^^ T2) [min(At 1 , r 2 ) - max(r 3 , 0)], (92) 

g? = (-iy +1 ^^[n + zr(r i )5 ll ], 1 = 0,1,2, < = 1,2,3,4, (93) 

where 

, , f 1 for x = , A 
ZT ^ = { for ' ^ 

and the definition 0(0) = \ is used again, we obtain finally the following closed-form ex- 
pressions for the representative BR geometric factors (6)-(8) with spherical space regions: 
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A ™ u) = t^rV E a^ ) E^°j i 4(o;i.i.o.';«i.«2.^^m(fl), (95) 

m 

4f } = T^I^V E 4- B) E # ji 4 (0; 1,1,0, 2; ifc, i? 2 , r„ R)Y 2m (R), (96) 
Bgp = (2 t)3 R.R^ E ^ JUM^); 1. 1, zr(r0-l, 1; fr, ifr, Ti, fl]yio(fl). (97) 
3. Numerical results and concluding remarks 

The convergence properties of the Fourier-Bessel expansions were examined in numerical 
calculations of some representative examples of BR geometric coefficients with spherical 
space regions, using the formula (67) with the Fourier-Bessel coefficients given by equations 
(68)-(73) and an expansion radius r ex = R + R 2 + max(T + At 2 , 0). 

For the BR geometric factors to have an appreciable magnitude, it is obvious that, 
in a system of units where the speed of light c = 1, the separation in space and time 
of the space-time regions must not be much greater than the dimensions of the space-time 
regions themselves. The calculations were done using an unspecified unit of length; choosing 
the millimetre as the unit of length, for example, the unit of time equals approximately 
3.33 x 10~ 12 s in a system of units with c = 1. This illustrates the fact that, on a realistic 
laboratory scale, the time intervals corresponding to even relatively large distances are very 
short, but we leave aside the question of how field measurements occupying and/or separated 
by such short time intervals can be realized. 

The numerical values of the representative BR geometric factors (6)-(8) for various 
spherical space-time regions with dimensions of the order of unity and similar or smaller 
space-time separations, in a c = 1 system of units, are collected in table 1. As the quantities 
that are required in the field commutation relations are the values of the differences Gjj' ir> — 
C^ l,l \ while some BR spring constants require the value of the sum C^' 1 ^ + C^' l \ the 
"reverse" geometric factor C^ 1 ' 1 ^ was calculated together with a given geometric factor C^ ,lV) . 
This was done by interchanging the radii R\ and i? 2 of the spherical space regions and 
changing their relative displacement R to —R (i.e., changing the polar angles 9 and of 
R to 7i — 8 and + 7r), together with interchanging the time intervals Ati and At 2 and 
changing their separation T to — T. BR geometric factors with fully coinciding space- 
time regions were calculated also as they are required for some of the BR spring constants. 
When the sphere separation R = 0, only the monopole (/ = 0) term contributes in the 
expansion (67) because ji(q n R) — > for / > as R — > 0, and thus the BR factors i^ 11 ' 
and -B^ 11 ^ with spherical space regions vanish when R = 0. The Fourier-Bessel expansions 
of the BR geometric factors, except those of the geometric factor i^ 11 ' for spherical space 
regions whose centres coincide, converge rapidly as less than a hundred terms were required 
for the four-digit accuracy with which the geometric factors are printed in table 1. 

The slow convergence of expansion (67) in the case of spherical space regions with co- 

inciding centres is due to the presence of a g-independent term in the quantity U 00 (q) of 
equation (68). However, the contribution of this term to the BR geometric factor A^ 11 ' has 
then a simple form and can be summed easily. With sphere separation R — 0, it follows from 
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equations (67) and (68) that the g-independent term contributes to the BR factor A^ 11 ' the 
quantity 

7ri?i i? 2 r ex Vr ex / Vr cx / 

where g^, defined by equation (90), is the q- independent term in question. The series in 
(98) is a Fourier-Bessel representation of of the simplest of the integrals evaluated in section 
3.2: 



— Ehi—RAjJ—Ra) = r miQM^q) dq 
r eX) ^ Vr ex / Vr cx / Jo 



ji 4 (0; 1,1, 0,0; i^ifc, 0,0) = ^. (99) 

Here, the result (82) is simplified using R < and _R>, which are respectively the lesser and 
the greater of the radii R± and R2, and the parameter r ex should be such that r ex > R > , 
which condition is guaranteed by that of equation (67b). When the result (99) is used in 
the expansion (67) in the cases of sphere separation R = 0, the rate of convergence improves 
dramatically and becomes similar to that of the geometric factors for space regions with 
noncoinciding centres. 

The accuracy of the numerical results of table 1, obtained using the Fourier-Bessel ex- 
pansions, was checked using the closed- form expressions (95)-(97). As no calculated values 
of the BR geometric factors could be found in the literature, this is the only check of the 
correctness and accuracy of our results. Admittedly, this check is not a fully independent 
one, as both the closed-form expressions and Fourier-Bessel expansions are obtained using 
Fourier-transform methods. However, the calculations reported here use only proven ana- 
lytical methods, of which Fourier integrals and Fourier-Bessel expansions are a part, and 
thus our results would be invalidated only if the same algebraic errors were made in the 
derivation of the closed-form and Fourier-Bessel expressions. In this connection, we note 
that the analytical expressions (48) and (49) for the double time averages, used in both the 
Fourier-integral and Fourier-Bessel formulations, were checked by doing the two-dimensional 
integrations involved numerically. 

A rather interesting result of the calculations is that the geometric factors A^ with fully 
coinciding space-time regions turn out to have negative values. The closed-form expression 
(95) for the geometric factor with radii R\ = R2 = Ro, time intervals Ati = At 2 = At 
and separations R = T = simplifies toi 

= -^( 4 + «)(2 - «0 2 ©(2 - «) - (100) 



8 Using this simple expression, it can be shown that, contrary to a conclusion of [H], the BR average 
'self-force' on the field measurement's test body approximates correctly the average self-force that 
obtains when the duration of the momentum measurements on a test body of sufficiently great 
mass is sufficiently short ||. 
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where k = Ato/Ro. For a fixed value of i?o, this function of the ratio k increases mono- 
tonically from — oo when k — > to the value of zero for k — > 0. For k > 2, the geometric 
factor iffl reduces to the value —1/RqK, and so the BR average 'self- force' |I| on the field- 
measurement's test body of charge density p\ is then pfVjAt D^A^' = —pjV^D^/Rl, 
which is simply the electrostatic force of attraction between the test and neutralizing bod- 
ies when their centres are displaced by a distance \Djp\ <C i?o- The negativity of the BR 
geometric factor A^ means that the spring constant fa = p\V?TiA™ of the spring that 
is used in a BR field measurement involving a space-time region Vi,Ti to compensate the 
test body's average 'self-force' has to be negative. While it is certainly possible to envisage 
a spring mechanism that would provide a force proportional to and in the direction of a 
test body's displacement,! we note that Bohr and Rosenfeld did not consider it necessary 
to make a comment on this rather unusual specification that their measurement procedure 
would place on the spring mechanism — but one can now only speculate whether Bohr and 
Rosenfeld were in fact aware of this consequence of their analysis. In any case, despite 
its inherent instability, a spring mechanism with negative spring constant should present no 
difficulty of principle for a BR measurement procedure because a BR spring, together with 
the test body to which it is attached, is supposed to be released only for the exact duration 
of the field measurement, and the spring force is designed so that its effect is compensated 
by the test body's 'self-force'. 

We conclude that a well-controlled method for the computation of the BR geometric 
factors was developed using Fourier-Bessel expansions. The efficiency and accuracy of the 
method were tested numerically in the case of spherical space regions when a BR geometric 
factor can be represented by a Fourier-Bessel series with terms expressed entirely in terms 
of elementary functions, and, using the computing system Mathematica, it is possible to 
obtain the factor in terms of manageable closed-form expressions. The space-time-averaged 
electromagnetic- field commutators, as well as the formal expressions and "gedanken" ex- 
perimental procedures of the famous Bohr-Rosenfeld analysis of the measurability of the 
electromagnetic field now can be translated easily and accurately into concrete numbers. 
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9 A spring system with a negative spring constant can be constructed as follows. Two elastically 
compressible rods, each of spring constant k and length I + d when not stressed, are aligned "head 
to tail" along a common axis and joined via a movable joint, while their outer ends are fastened 
using similar joints to rigid supports so that this spring system is compressed to a total length 21. 
It can be shown easily that a body attached to the middle joint and moved a small distance x, 
x <C I, x <C d, from the system's axis and in a direction perpendicular to it, will experience a force 
F = {2kd/l)x = —kx, which is proportional to and acting in the direction of the displacement x, 
i.e., the spring constant k = —2kd/l of such a system is negative. 

10 A hint that Bohr and Rosenfeld were aware of the possibility of the geometric factor A£x being 
negative is given by their careful writing of its square root as l-Ajj^l 1 / 2 . 
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the author to realize that the formal expressions for some BR geometric factors imply a 
suitable regularization. 



TABLE I. Representative BR geometric factors Cfj' for space-time regions I and II specified 
by space spheres of radii R\ and R2, and time intervals At± and At 2 , respectively, with the centre 
of the second sphere displaced from that of the first one by a vector of spherical coordinates R, 
9, 4>, and the beginning of the second time interval separated from that of the first one by a time 
interval T; units such that the speed of light c = 1 are used. 





Ri 


R2 


R 


6 




At 1 


At 2 


T 




1(1,11) 

■flu 


l 


1 









1 


1 





-1.625 x 10+° 


I(I,II) 


10 


10 









1 


1 





-2.850 x 10- 3 


I(I,II) 


1 


1 









1 


2 


0.5 


1.953 x 10- 1 


l(I,H)a 


1 


1 









2 


1 


-0.5 


-5.664 x 10- 1 


j(i,n) 

-rt-XX 


1 


1 




h 


¥ 


1 


1 


0.5 


-6.407 x 10~ 2 


l(I,H)a 
si-xx 


I 


1 




h 


h 


1 


1 


-0.5 


-4.530 x 10" 1 


I(I,II) 

■t^xy 


I 


1 




h 


h 


1 


1 


0.5 
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-8.914 x 10~ 2 
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3.493 x 10~ 3 
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1 
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-3.884 x 10" 4 
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o(I,H)a 
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1 
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4.126 x 10 -3 



a The "reverse" geometric factor Cjj ' ' of the preceding entry. 
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